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2^ ■ Abstract 



Many authors - beginning with Bekenstein - have suggested that the energy 



O 
O 

(3 ■ levels En of a quantized isolated Schwarzschild black hole have the form 

^ ■ En = (Jy/nEp, n = 1, 2, . . . , 0" = 0(1), with degeneracies g"'. 

. In the present paper properties of a system with such a spectrum, considered 

^^H, as a quantum canonical ensemble, are discussed: 

5^ ' Its canonical partition function Z{g,f] = 1/kBT), defined as a series for 

g < 1, obeys the 1-dimensional heat equation. It may be extended to values 
g > 1 hj means of an integral representation which reveals a cut of Z{g,(3) 
^ . in the complex g-plane from g = 1 to g oo. Approaching the cut from 

above yields a real and an imaginary part of Z. Very surprisingly, it is the 
(explicitly known) imaginary part which gives the expected thermodynamical 
properties of Schwarzschild black holes: 

Identifying the internal energy U with the rest energy Mc^ requires (3 to have 
the value (in natural units) 

(3 = 2M(ln g/a^) [1 + 0(1/ M^)] 

[A'Ka'^ = lng gives Hawking's jSn) and yields the entropy 

S = [In g/iAna^)] A/ A + 0(ln A) , 

where A is the area of the horizon. 
-'^ E-Mail: kastrup@physik.rwth-aachen.de 
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1 Introduction 



As the (Hawking) temperature 

(Ep = ^c^h/G is Planck's energy) of the radiation emitted by a black hole[Q 
is proportional to Planck's constant h, i.e. a quantum effect, it was clear from 
the beginning that its deeper understanding would require a quantum theory 
of the gravitational properties of black holes. Despite the possible lack of a 
convincing general quantum theory of gravity many attempts have been made 
to identify the quantum energy levels En of an isolated Schwarzschild black 
hole. Bekenstein[^ was the first to use Bohr-Sommerfeld type quantisation 
arguments and suggested a spectrum 

= av^Ep ,n = 1,2,..., (2) 

where a is a (model dependent) dimensionless constant of order 1. 

Since then a number of authors[3-23] have given different arguments for a 

quantum black hole spectrum of the type (2). 

I myself discussed in ref. [19] how such a spectrum may be understood in the 
framework of a stringent canonical quantisation of the purely gravitational 
Schwarzschild spherically symmetric system P^||25||. 

In the following I shall take the relation (2) for granted and ask what its 
implications for the thermodynamics of the system are if this is viewed as 
a canonical ensemble. Again I shall not enter into a discussion of possible 
conceptual problems of such a statistical framework in the context of black 
holes[26-29,31] and shall make a few comments after the results have been 
presented. 

As to the present status of the thermodynamics of black holes see the reviews 
by Wald[p^, |30[| , Brout et al.[^] and Sorkin[|32|. The present situation as 



to the quantum states of black holes in the framework of string theories 
has recently been reviewed by Horowitz . (Strings can have a spectrum 
171^ oc n, too!) 
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2 Canonical partition function 



Following the recent discussion by Bekenstein and Mukhanov|jT5[] (and that 
by the same and other authors before) I assume the degeneries dn of the 
levels (2) to be (those authors actually take dn = 2"~^ which is, however, 
not essential for the arguments below). It is convenient in the following to 
define t = \ng. 

The canonical partition function of the system is 

Z{t,x) = J2e^''^e~^^^ = l + e^Z{t,x) , t = Ing , x = (3aEp , (3) 

n=0 

where Z is the partition function corresponding to the assumptions of Beken- 
stein and Mukhanov. 

The series (3) is obviously convergent for t < {\g\ < 1) and converges 
for t = {g = l),x > 0, according to the Maclaurin-Cauchy integral 
criterium||34||. For t > the series is divergent, but the function Z{t,x) 
can nevertheless be defined by continuation (see below). 
The series (3) obeys the heat equation 

dtZ = dlZ , (4) 

with the boundary values 

oo . 

Z{t ^-oo,x) = 1 , Z{t = 0,x) = Y.e~^^^ = (p{x) , (5) 

n=0 

Z{t,x -^oo) = 1 , Z{tj^O,x = 0) = ^-^ = r]{t) . (6) 

1 — e'' 

For the series 0(x) we have the lower and upper bounds 

" due-^'' = 2(- + i^)e-^ < 0(x) -1< T rfz/e"^^ = ^ . (7) 
X Jo x^ 

Observing the behaviour of e~^^^ between n = and n = (/e + 1)^, {k + 
1)^ — /c^ = 2/c + 1 one can sharpen these bounds] 



(E2te-*-)-e- = (8) 



< 0(x)-l<^2fce-^^= _^ 
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As to the physics of the system one is interested in properties of the partition 
function for t > 0. At a first sight one might consider to solve the heat 
equation for t > 0,x > with the known initial values (5) and (6) in a 
standard manner! 
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Zit,x) = dy[K{t,x-y)-K{t,x + y)](t){y)+ j dTK{t - t,x)7]{t) , (9) 
where K is the "heat kernel" 

K{t,x) -- 
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and K, essentially, its x-derivative 

k{t,x 



^e-^V(4t) 



2V7rt3 



-2d,K{t,x) . 



(10) 



(11) 



This approach does not appear to work, however, because the functions (p 
and 1] are not "decent" enough p^, 37]: rj becomes singular as 1/t for t — »• 
and behaves like for x — ^ 0. The latter property can be inferred from 
the inequahties (7) and (8). 

An extension of the function Z{t,x) - actually it was Z - defined by the 
series (3) into the complex gf— plane was discussed 100 years ago by the 



mathematician Lerch 38 : Using the relation [39 



nx^ 



nr Jo 



(12) 



\x\ 



2J^Jo 



dr 



-x^/(4r) — riT 



drK^T, x)e 



-nr 



converts the series (3) into a geometrical one which can be summed under 
the integral for t < 0, a; > 0: 



Z{t = \ng, x) 



71 Jo 
X 



l-e(^-lK) 



20F Jo r3/2 

X , e' 

du— 



~ dr -xV(4r)_ 1 



2^Ji 



-X? I (4 Inw 



1-e 
) 1 



In^/^ 



u 



u- g 



(13) 
(14) 

(15) 
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Notice that the relation (14) may also be written as 

Z{t,x)= K{T,x)r]{t-T) = dfK{t-f,x)ri{f) , (16) 



where t = t — t. 
Observing that 

y/TC Jo 

we get 

Z-l = e'Z{t,x) = e*^ / rft;e-^^/4 ^- (17) 



e* / dTK(T,x)e~^r](t-T) . 
Jo 



The expressions (13) etc. may also be obtained by inserting the representation 
(12) of exp{—^/nx) into the series (3) and Borel summing ||4CI|| it for t < 0. 
The integrals converge for all values of g l,t ^ (the ones representing 
Z converge even better than those for Z). For real g > l{t > 0) one has to 
take the principal value of the integrals. 

The integral representations (13)-(15) for Z(t,x) are solutions of the heat 
equation for all (even complex) t 7^ as can be seen immediately, e.g., by 
replacing the differentiation of 1/(1 — e''' ~ ^) in eq. (14) with respect to t 
by the negative one with respect to r and performing a partial integration 
afterwards, or from the eq. (16) directly, because K{t,x) is a solution of the 
heat equation. 

Notice that Z(t,x) is not a solution of the heat equation (4), only e^Z(t,x) 

is one. 

As 

Z{t,\x) = ^ dve-'"''/^ , , A>0, (18) 

Vvr Jo gtg-A /v 

we see that 

lim Z{t, Xx) = — ^ , lim Z{t, Xx) = 1 , (19) 
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in accordance with eqs. (6), but now for t > 0. 
Correspondingly we get for Z(t,x): 

lim Z{t, \x) = , lim Z{t, \x) = e~^^ . (20) 

The expressions (13)-(15) or (17) can be used to extend the function Z{t = 
\ng,x) or Z{t,x) into the complex g- or t-planes: 

According to eq. (15) Z{g,x) has a branch cut in the complex g-plane along 
the real axis from 1 to oo. The discontinuity of Z across the cut is given 
by|8|, P 

lim [Z{g + ie) - Z{g - ie)] = 2mk{t, x) , (21) 
e ^ 0+ 



and Z{g,x) is an analytic function of g except for this cut[^, ^ 



If one, therefore, approaches the cut along the real axis from above, the limit 

lim Z{g + ie,x), g> I , 
e ^ 0+ 

is no longer a real-valued function of g but has a nonvanishing imaginary part 
Zi{t, x). Standard procedures used in the field of dispersion relations[^] yield 
as the real part Zr{t,x) the principal value integral 

POO 

Zr{t, x) = p.v. / dr K{t, x)7]{t - t) , (22) 
Jo 

and for the imaginary part Zi the expression 

Z.(t,x)=vri^(t,x) = ^e-^V(4t) . (23) 

Obviously Z^ and Zi are solutions of the heat equation seperately. Zr{g,-) is 
the Hilbert transformer], Q of K{g, ■). 

We shall see that, strangely enough, it is this imaginary part of the partition 
function which gives exactly the thermodynamical properties expected for 
black holes! 

The formal reasons for this can be seen from the limits (19) which obviously 
are those of Z^: The limit of Zr{t, x) for x — > is negative for t > 0, whereas 
Zr is positive for large x\ Notice that the integrand in eq. (22) is negative for 
t > r > and that, therefore, Zr may be negative for small x. A negative 
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partition function can, however, hardly be interpreted thermodynamically 
where the logarithm has to be taken. On the other hand the imaginary part 
(23) is positive for alH > 0,x > and does not have the " sign-desease" of 



3 Thermodynamics 

In the thermodynamics of black holes one is especially interested in the be- 
haviour of the system for large /5 (low temperatures), because the inverse 
Hawking temperature = l/(A;BTff) is very large for macroscopic black 
holes (see eq. (1)). 

According to cq. (20) the real part Zr of the partition function Zit^x) be- 
haves like exp(— x) for large x and therefore the associated internal energy 
U = —din Zr/d(3 is just Ep, i.e. the lowest possible energy level, as one would 
expect naively from the conventional paradigms of statistical mechanics. 
The main reason for not using the real part Z^., however, is its property to 
become negative for small x ii t > as discussed above in connection with 
eq. (22). 

The situation becomes surprisingly interesting and unconventional if we take 
the imaginary part Zj, eq. (23), as the partition function for calculating the 
thermodynamical properties of the system: 
We first calculate the internal energy 

[, = £ = -^ = -^=M.^ . 

dp Zi 

which we identify with the total rest energy Mc^ of the black hole and get 

U^^P-^ = Mc' . (24) 
Solving this equation for (3 and discarding the negative root gives 

which for Mc^ 3> Ep leads to 



2tM(? a^E"^ 



/3=-^(l + ;^TT7£l)- (26) 
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For = t/(47r) we get Hawking's Ph, plus a small correction of order 

iEl/iM'c'))f3H\ 

Furthermore, for the average of the energies squared we have 

^ = ^'^'^^/^^ = = - YtyEl , (27) 

yielding the mean square fluctuations 

2t /5 



We see that for t > the right hand side is negative, as expected, corre- 
sponding to a negative specific heat. Considering the fact that the left hand 
side of eq. (28) is the difference of two very large numbers, both of the order 
M^c^, the mean square fluctuations are actually very small and essentially 
universal, because the right hand side depends on M itself only through the 
term which is negligible for very large (3\ 

If we define the average level number N by (see eq. (3)) 

N-^/Z., (29) 

then we get 

iV = (^ - A) = W/{a'El) ^ {Mcr/{a'El) . (30) 

The last approximate equality follows from the fact that the fluctuations (28) 
are so small. 

Thus we have for large N 

X = aEpf3 = 2tVN = 2t—— . (31) 

(7 Ep 

Finally we come to the entropy 

S/kB^lnZi + (3U , (32) 
for which we get the exact result 

S/kB^^ + lnx-^\nt + \n{y/^/{2e)) , (33) 



or, using the relation (31) and ignoring terms of order 0(1), 

S/kB = iTV+^lnTV (34) 

where A = 47ri?| = IGttG^M^/ is the area of the horizon with Schwarzschild 
radius Rs and /p = Gh/c" the Planck length squared. 

For cr^ = t/(47r) the leading term of S/ks has the Bckcnstein-Hawking value 
74/(4Z|,), for (J = 1/2, as discussed in ref. [19], it is slightly smaller if i = ln2 
etc.. 

Notice that the factor exp(— a;^/4t) - typical for the heat equation - is the 
decisive one for providing the essential features of the thermodynamical prop- 
erties just discussed. 

Up to now we have been dealing mainly with the large (3- (low temperature)- 
behaviour of the system. There arc in addition some features for small j3 
(high temperature) worth mentioning. Whether they have any physical sig- 
nificance - e.g. for the big bang era - remains to be seen. 
The internal energy U in eq. (24) vanishes for x = xi = \p2A or ksTx = 
aEp/V2i. W becomes negative (eq. (27)): E^i = -a'^El/tl The entropy 
(33) at this temperature takes the value Si/kB — (1/2) ln[7r/ (2ei^)] which for 
t = In 2 is equal to 0.092 ■ ■ ■! 

According to eq. (27) the quantity E"^ vanishes for x = X2 = ^/Qi,kBT2 = 
aEp/^M < fcsTi, with U = U2 = 2aEp/\f&i,S2/kB = (1/2) In[3e7r/(2t2)] > 
SjkB. 

Then there is the value x — Xo{t) for which the entropy (33) vanishes. The 
resulting equation 

cannot be solved for Xo{t) explicitly, but the inequality Xo{t) < 2et^l'^ j \pK 
follows immediately. 

Strange things seem to happen at temperatures of the order of the Planck 
energy! 
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4 Remarks 



The most surprising feature of the "canonical" quantum statistical mechanics 
of the level spectrum (2) is that the expected thermodynamical properties 
of a Schwarzschild black hole are associated with the imaginary part of the 
partition function for real g > 1\ 

Superficially this appears to question the canonical approach to the thermo- 
dynamics of black holes as possibly inappropriate. On the other hand the 
thermodynamical properties associated with Zi are too intriguing and too 
interesting in order to dismiss them. The formal reasons for preferring Zi 
compared to Z^ as the partition function relevant for the thermodynamics of 
the system have already been stressed in connection with the eqs. (22) and 
(23). The "physical" thermodynamical consequences fully justify the more 
formal conclusions! 

There is in addition an interesting heuristic consistency argument why the 
imaginary part (23) of the partition function is the physically relevant one: 
Suppose the sum (3) does not extend up to oo but up to a very large number 
A^. Then we get 



instead of eq. (14) and the integrand is not singular anymore for t = t (the 
fraction resulting from the finite geometrical series just has the value for 
T = t) and there is no cut! If we now exploit eq. (31) and replace the x^ 
under the integral sign by At'^N we get 



Jo 



oo 



k{r) , (37) 



Kr)=t-r--^ k{r) = -^ 1 _ ,(t - r) 



Evaluating the integral for large N by means of a saddle point approximation [|3 
yields a saddle point for r = t and the result 

Z„ « |ViVe-«A' = |^e-V(4i) , (38) 
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where again the relation = At'^N has been used. 

The approximation (38) is not equal to the function (23) and not a solu- 
tion of the heat equation - there is mainly one additional factor x/^/t due 
to the approximations involved - but it contains the most essential factor 
exp(— a:;^/(4t)) which is so important for the qualitative structure of the ther- 
modynamics. 

Furthermore, the smallness of the thermal fluctuations (28) show that the 
thermal interactions of the black hole with the heat bath are small and, 
therefore, a "canonical" statistical treatment seems plausible and may not 
be too far off a microcanonical one. 

In any case, the above "canonical" results have to be interpreted as require- 
ments on the properties of a heat bath if it is to be in thermal equilibrium 
with the black hole. 

If the whole approach discussed above is not unsound then the level spectrum 
(2) has to be taken seriously and a convincing justification of its validity is 
desirable fii I . 

A final remark: The thermodynamics discussed would be quite different if 
we would not interpret g = exp(t) as a fixed number but as the tempera- 
ture dependent fugacity g = z = exp(/i/3), /x: chemical potential, of a grand 
canonical ensemble [E^ . 



I thank C. Gutsfeld, G. Roepstorff, F. Schramm, T. Strobl and H. Wissowski 
for helpful discussions and critical remarks. Especially the numerical tests 
of C. Gutsfeld and F. Schramm together with the joint discussions including 
T. Strobl prevented me to expect too much from the real part Zr(t,x), eq. 
(22), and stimulated me to take the imaginary part Zi more seriously. 
Finally I thank my wife Dorothea for her support, her understanding and 
her patience while this paper was being prepared. 



Note added: 

After this paper was submitted as an e-print M. Perry kindly drew my at- 
tention to refs. pB| , where the imaginary part of the partition function is 
related to the metastable states of the system. Especially the last paper [47] 
is of considerable interest in the context of the present article. 
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